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INTRODUCTION 


The  classification  of  sonar  contacts  as  submarine  or 
non-submarine  is  a primary  concern  to  Naval  forces  assigned  an 
ASW  mission.  Incorrect  classification  results  in  undue 
expenditure  of  resources  to  investigate  or  conduct  an  attack  on 
a non- submarine,  or  in  permitting  a true  submarine  to  proceed 
unchallenged.  Though  the  classification  problem  at  short  sonar 
ranges  has  not  been  solved,  the  effects  and  sources  giving  rise 
to  short  range  false  contacts  have  received  considerable 
attention  and,  in  most  cases,  can  be  described.  Sufficient 
experience  in  the  use  of  Bottom  Bounce  (BB)  and  Convergence  Zone 
(CZ)  modes  of  propagation  has  not  been  accumulated  to  adequately 
evaluate  the  character  and  extent  of  the  classification  problem 

with  long  range  sonars  using  these  propagation  paths. 

Ti  s 

' study-^was  initiated  at  TRACOR,  Inc.,  in  March  1969 

under  Contract  No.  N00024-69-C-1255  to^evaluates  the  character 
and  extent  of  the  long  range  classification  problem.  It  was 
decided  to  investigate  first  the  BB  mode  of  propagation, 
reserving  the  CZ  mode  for  future  work.  An  investigation  by 
TRACOR^has  shown  that  schools  of  fish  are  probable  sources  of 
false  contacts  for  sonars  operating  with  a CW  pulse  in  the  BB 
mode.^. 

v ’ — — - j 

For  the  purpose  of  that  investigation  a false  target 
was  considered  to  be  any  non- submarine  producing  an  echo  whose 
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characteristics  caused  it  to  be  classified  as  a submarine  echo 
on  the  basis  of  echo  strength  and  duration.  For  a school  of  fish 
to  be  mistaken  for  a submarine  several  conditions  must  exist. 
First,  the  target  strength  of  the  school  must  be  comparable  to 
that  of  a submarine.  Second,  the  extension  of  the  school  in  the 
direction  of  sound  propagation  should  not  be  greater  than  the 
length  of  a submarine  so  that  the  echo  will  not  be  longer  than 
that  expected  from  a legitimate  target.  Finally,  since  the 
azimuthal  beam  width  at  ranges  of  interest  in  bottom  bounce 
operations  is  much  greater  than  the  length  of  a submarine, 
detection  on  more  than  one  beam  is  highly  unlikely.  Hence,  the 
school  of  fish  should  subtend  an  azimuthal  angle  no  greater  than 
the  azimuthal  beam  width  of  the  sonar. ■ 

Ray  theory  was  used  to  determine  the  volume  which  could 
be  occupied  by  the  fish  school  consistent  with  the  above 
constraints.  The  target  strength  for  the  school  was  then 
calculated  by  multiplying  the  target  strength  of  a single  fish 
by  the  number  of  fish  in  the  volume.  This  procedure,  of  course, 
neglects  all  multiple  scattering  effects.  The  target  strength  of 
a single  fish  was  taken  from  measurements  made  by  Cushing,  et.  al. 
It  was  found  that  schools  of  fish  could  exist  which  have  a target 
strength  comparable  to  that  of  a submarine. 

Now  it  is  known  that  some  schools  of  fish  are  so  dense 
that  the  fish  are  in  actual  physical  contact.  One  would  expect 


for  such  a school  that  multiple  scattering  effects  could  not  be 
ignored.  If  this  is  the  case,  then  classification  characteristics 
such  as  (1)  echo  duration,  (2)  spectrum,  (3)  sharpness  of  onset 
and  trailing  edge,  (4)  highlight  structure,  (5)  doppler,  (6)  target 
strength  and  (7)  ping-to-ping  consistency  and  persistence  would 
differ  from  those  calculated  on  the  basis  of  neglecting  multiple 
scattering.  Thus,  since  this  initial  investigation  showed  that 
fish  schools  posed  a problem  to  long  range  classification,  it 
was  decided  that  the  only  way  to  adequately  evaluate  the 
character  and  extent  of  this  problem  was  to  formulate  a realistic 
model  of  the  multiple  scattering  of  a pulse  from  a school  of  fish. 
This  preliminary  report  takes  the  first  step  in  developing  such 
a model. 

In  Part  I,  the  general  theory  of  the  multiple  scattering 
of  an  acoustic  pulse  from  a random  collection  of  point  scatterers 
is  developed  in  some  detail.  Only  the  first  order  statistics  of 
the  scattered  signal  will  be  considered  in  this  report.  The 
second  order  statistics,  which  are  related  to  the  scattering  cross 
section  of  the  fish  school,  will  have  to  be  reserved  for  future 
work. 

In  Part  II,  the  general  theory  developed  in  Part  I is 
applied  to  the  case  where  the  random  collection  of  point  scatterers 
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represents  a dense  school  of  fish.  The  scattering  properties 

3 

used  for  a single  fish  are  those  reported  by  D.  E.  Weston  . 

Thus,  it  is  assumed  that  the  scattering  is  due  entirely  to 
the  air  bladder  of  the  fish  and  that  the  fish  can  be  represented 
as  an  isotropic  point  scatterer.  This  will  be  approximately  the 
case  for  a fish  of  length  L if  the  frequency  of  the  sound  wave  is 
less  than  cq/L,  where  cQ  is  the  speed  o£  sound  in  water.  For 
frequencies  greater  than  cq/L  scattering  from  the  fish  tissue 
becomes  increasingly  important  and  the  simple  theory  of  isotropic 
scattering  from  an  air  bladder  does  not  apply.  The  theory  of 
scattering  from  bladder  fish  takes  into  account  dispersion 
and  absorption  in  the  scattering  volume. 

Two  cases  are  considered  for  the  geometry  of  the 
scattering  region  simply  bece-  <■»:  they  are  the  most  amenable 
to  solution:  (1)  the  fish  are  contained  within  a spherical 
region;  (2)  the  fish  are  contained  within  a layer  with  infinite 
plane  boundaries.  In  the  first  case,  spherical  schools  of  fish 
have  been  observed  in  the  open  sea,  but  the  diameter  of  the 
schools  and  the  distribution  of  fish  in  them  were  not  reported. 

In  the  second  case,  while  a layer  with  plane  boundaries  is 
amenable  to  solution,  it  does  not  meet  all  the  requirements 
imposed  on  fish  schools  which  were  discussed  at  the  beginning 
of  this  section.  It  does  not  meet  the  requirement  that  the 
school  be  contained  entirely  within  a single  beam  since  the 


plane  boundaries  of  the  layer  extend  to  infinity.  However,  this 
would  serve  as  a good  model  for  the  Deep  Scattering  Layer  (DSL) 
since  there  is  a growing  amount  of  evidence  to  indicate  that  the 
DSL  is  composed  of  bladder  fish. 

. At  this  time  it  is  anticipated  that  the  first  order 
statistics  for  the  scattered  field  will  be  numerically  evaluated 
for  the  spherical  school.  However,  no  further  work  is  planned 
on  the  model  for  the  DSL. 
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Part  I.  MATHEMATICAL  THEORY  OF  MULTIPLE  SCATTERING 
A.  INTRODUCTION 

The  problem  of  scattering  of  an  incident  wave  by  a 
single  obstacle  has  been  considered  rather  thoroughly  in  the 
literature,  beginning  with  Rayleigh's  work  in  fluids,  and 
continuing  to  the  present  day  with  various  quantum-mechanical, 
electromagnetic,  and  elastic  cases.  The  literature  on  multiple 
scattering,  on  the  other  hand  is  not  so  extensive.  Regular 
arrays  of  scatterers  were  treated  by  Huygen's  principle  or  various 
perturbation  schemes,  mainly  with  an  eye  to  obtaining  their 
"strong"  filtering  properties  arising  due  to  periodicity. 

Out  of  this  work  came  X-ray  diffraction  theory  and  the  band 
theory  of  solids.  Until  1945  multiple  scattering  from  a random 
distribution  of  scatterers  was  restricted  to  the  scattering  of 
particles  based  on  the  Boltzmann  integro  differential  equation 
describing  transport  processes.  This  formulation  is  merely  the 
expression  of  conservation  of  particles  in  phase  space;  hence 
the  treatment  is  classical,  with  no  account  taken  of  the 
quantum-mechanical  wave  nature  of  the  particles  or  photons. 

Such  a theory  would  be  expected  to  be  valid  only  it  the  wavelength 
of  the  particles  is  much  smaller  than  the  average  distance  of 
separation  between  the  scatterers. 

There  are  also  a large  number  of  problems  of  multiple 
scattering  in  which  the  wavelength  is  comparable  to  the  average 
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scatterer  separation;  some  examples  of  this  latter  type  of 
problem  are  acoustic  wave  propagation  in  bubbly  water,  scattering 
from  a school  of  fish  and  the  scattering  of  electrons  or  x-rays 
by  the  nuclei  of  liquids  or  amorphous  solids.  Any  treatment  of 
these  problems  must  include  the  reflection,  refraction,  and 
interference  phenomena  that  are  characteristic  of  wave  problems; 
hence  it  must  be  based  on  the  wave  equation,  rather  than  on  the 
simple  conservation  statement  leading  to  the  Boltzmann  equation. 

The  first  systematic  treatment  of  the  problem  of  the 
multiple  scattering  of  waves  from  a random  distribution  of 
scatterers  was  published  by  Foldy^  in  1945.  Foldy's  unique 
contribution  was  the  introduction  of  the  concept  of  "configurational" 
averaging  of  relevant  physical  quantities  by  defining  a joint 
probability  distribution  for  the  occurrence  of  a particular 
scatterer  configuration.  By  averaging  the  equations  of  multiple 
scattering  over  the  statistical  ensemble  of  scatterer 
configurations,  Foldy  was  able  to  derive  integral  equations 
governing  these  configurational  averages. 

Foldy's  work  dealt  only  with  a random  collection  of 
isotropic  point  scatterers.  This  was  later  generalized  by  Lax"* 
to  treat  the  multiple  scattering  of  quantum  mechanical  waves 
by  point  scatterers  having  quite  general  scattering  characteristics, 


£ 

and  by  Waterman  and  Truell  to  treat  scattering  regions  having 
non-vanishing  dimensions.  A unified  picture  of  the  propagation 
of  the  coherent  and  incoherent  radiation  associated  with  multiple 
scattering  is  given  by  McCloskey^.  Surveys  of  the  entire  field 

of  the  multiple  scattering  of  waves  have  been  published  by 
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Twersky  , and  Burke  and  Twersky  . 

To  the  best  of  the  author's  knowledge,  the  theory  of 
multiple  scattering  from  a random  collection  of  scatterers  which 
has  appeared  in  the  literature  to  date  has  only  been  concerned 
with  the  scattering  of  an  incident  monochromatic  wave.  In  this 
report  an  attempt  is  made  to  extend  the  theory  to  include  the 
scattering  of  a pulse  (time  limited  signal).  Only  the  first 
order  statistics  of  the  scattered  signal  will  be  considered  here. 
The  second  order  statistics  of  the  scattered  signal  will  be 
considered  in  a later  report. 

In  Part  I a detailed  description  of  the  theory  of 
multiple  scattering  of  a pulse  from  a random  collection  of 
volume  scatterers  will  be  presented.  Results  found  in 
References  4-7,  mentioned  above,  will  be  freely  used  without 
any  additional  references  to  these  papers  being  made. 

B.  STATISTICAL  PRELIMINARIES 

Assume  we  have  a collection  of  N point  scatterers  and 
that  the  iL  scatterer  can  be  characterized  by  its  position 
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vector  r^  = (x^,y^,z^)  and  a scattering  parameter  a^.  Here 

x. ,y.  and  z.  are  the  Cartesian  coordinates  of  the  1 scatterer. 
i,Jx  i 

The  only  application  that  will  be  considered  is  the  scattering 
from  bladder  fish.  For  this  case  it  is  assumed  that  the  scattering 
is  due  entirely  to  the  swim  bladder  of  the  fish  and  that  the 
actual  elongated  bladder  can  be  replaced  by  an  equivalent  spherical 
air  bladder.  Consequently,  the  scattering  parameter  a^  will  be 
taken  to  be  the  radius  of  the  equivalent  spherical  bladder.  There 
is  no  difficulty  in  generalizing  the  theory  if  more  than  one 
scattering  parameter  must  be  associated  with  each  scatterer. 

We  shall  say  that  we  have  a particular  configuration 

-> 

of  scatterers  if  the  position  vectors  r^  and  the  scattering 
parameters  a^  are  specified  for  each  of  the  N scatterers.  This 
configuration  may  be  regarded  as  one  state  in  an  ensemble  and 
an  average  over  all  states  may  be  taken.  In  order  to  simplify 
the  notation  in  what  follows,  a four  dimensional  vector  q^  will 
be  introduced  and  defined  to  have  as  its  components  q^  = 

(ri,ai^  = (xi»yi»zi»ai)*  Also,  will  represent  a four 
dimensional  volume  element  and  will  be  defined  as  dcj^  = 
d?. da.  = dx. dy. dz . da . . Thus,  we  can  now  say  that  we  have  a 
particular  configuration  of  scatterers  if  the  vectors 
q^  (1<  i<N)  are  specified.  This  four  dimensional  space  will  be 
called  phase  space. 


!! 
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To  further  specify  this  ensemble,  we  introduce  a 
joint  probability  density  function  p(q^>  • • • >q^)  so  that 

P(qi»  • • • • *dqN 

represents  the  probability  of  finding  the  scatterers  in  a 
configuration  in  which  the  first  scatterer  lies  in  an  element 
of  volume  dr^  about  the  point  r^  and  has  a scattering  parameter 
lying  between  a-^  and  a^+da-^,  while  at  the  same  time  the  second 
scatterer  lies  in  an  element  of  volume  d^  about  the  point 
and  has  a scattering  parameter  lying  between  &2  anc*  a2+da2»  etc. 
The  integral  of  the  joint  probability  distribution  over  all 
configurations  of  scatterers  is  normalized  to  unity;  i.e., 

J...S  p(4i> ■ • • ,?N>  tffr • -dqN  ‘ 1 • (i) 

We  can  also  introduce  the  number  density  n(cfp  . . . ,c^  ) of 
scatterers  by  multiplying  through  Eq.  1 by  N and  letting 

n(q^ » • • • ) = N p(q^>  • • • >q^)  • ^ ) 

The  integral  of  n(q^, . . . ,q^ ) over  all  of  phase  space  gives 
the  total  number  of  scatterers  present;  i.e., 

n(q^> •••>%)  d^l*-*d^j  = N * (3) 

The  probability  of  a configuration  having  the  iL 
scatterer  occupying  the  volume  element  dq^  regardless  of  the 
location  and  value  of  the  scattering  parameter  of  all  the  other 
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scatterers,  may  be  obtained  by  integrating  over  all  but  the 
i coordinate  as  follows: 

p(q^)  =/.../  p Cq ^ ) dqi*  • -^i-i^i+i*  • *dqN  * 

In  terms  of  the  number  density  we  can  write 


N p(q^)  = n(qL)  . (5) 

It  is  also  convenient  to  introduce  the  notion  of 
conditional  probabilities.  The  conditional  probability, 

p(q^, . . . ,qt  _1,  q^+1, . . . ,cfN|  if.^)  , for  a configuration  having  the  iC 
scatterer  fixed  at  the  known  location  “r^  and  having  a scattering 
parameter  with  the  known  value  a^  is  defined  to  be 

P(ql’ • • • >qi-l>^i+l> ' * * iqi'  ^ ^ 

For  example,  if  i = 1,  Eqs.  4 and  6 become 
p(qi)  =/•••!  dq2---d?N  > (?) 


p(q2>  • • • > qj  ql>  =7^7 


p(q1} • • • ,qN)  • 


These  notions  can  be  generalized  to  the  case  where 

more  than  one  scatterer  is  held  fixed.  To  facilitate  writing 

let  us  take  the  case  where  scatterers  1 and  2 are  held  fixed 

th  th 

instead  of  the  more  general  case  where  the  1 and  j scatterers 
are  held  fixed.  First,  the  joint  probability  of  finding  the  first 
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scatterer  in  the  volume  element  dq^  and  the  second  scatterer  in 
the  volume  element  dc^  is  given  by 


P (q  X > ^2  ) J*  • • • J*  P » • • • j ^ ^q3 ' ' ' ^9  i 


(9) 


3 "N 

The  conditional  probability  for  a configuration  having 
the  first  scatterer  fixed  at  ^ and  having  a scattering  parameter 
with  a value  a-^  and  the  second  scatterer  fixed  at  r ^ and  having 
a scattering  parameter  with  a value  Is  defined  as 

p(qi»---»qN)  . 

p(q1}q2)  (10) 


p(q3>  • • • - 


Another  useful  relation  can  be  obtained  from  these 
results.  Combining  Eqs.  8 and  10  gives 

r7 

->  , P<~ 

I di .a« ) = 

N 


p(q3,  • • • ti,q2) 


p(?2,  • • . >3n|  q\)p  C?! ) 


P(q!,q2) 


(11) 


Moreover,  Eq.  8 for  the  case  N=2  yields 


1 (^2 1 ^1 ) ~ 


p(qi>q2) 


p(qx) 


(12) 


Solving  Eq.  12  for  pA-p^)  anc*  substituting  this  result  in 
Eq . 11  gives  the  desired  result 


p(q2»- • • »’qN|’?x)  = p(53»---7qNfqi»t2)  p(q2fqx)  . 


• • • 


(13) 


As  a final  illustration  of  the  use  of  number  densities, 
we  can  write  Eq  12  as 


p(q\,l2>  * p(qi>p(q2!  <£>  - — 


(14) 


Here  the  conditional  number  density  n ( q 2 1 ) has  been  introduced. 

The  conditional  number  density  n{q"2|q\)  is  related  to  the 
conditional  probability  density  pC^jq^)  in  the  same  way  that 
n(q^)  is  related  to  p(q^)  (cf.  Eq.  5);  i.e.,  by  the  relation 

n(t2|ti)  = P(?2|ti>  * ‘ (15> 

The  factor  (N-l)  appears  in  Eq.  15  simply  because  one  particle 
is  held  fixed. 

The  generalization  of  these  results  to  the  case  where 
an  arbitrary  number  of  scatterers  are  held  fixed  with  known 
values  of  their  scattering  parameter  is  immediate. 

Now  consider  a function  f (?,  t fq^, . . . »<^)  of  space 
coordinates  r,  time  t and  the  N vectors  q^,...,q^.  The 
configurational  (or  ensemble)  average  of  f (r,  t|  c^, . . . ) is 

defined  by 

(f(r,t^>  =J.  . .Jf(r,  tj^, . . . ,^N)  p(q^, . . . ,tN)  d^.  . .dqN  . (16) 
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The  partial  average  with  one  or  more  scatterers  held 
fixed  is  obtained  by  averaging  over  the  appropriate  conditional 

probability:  Thus,  the  average  of  f(?, t | qN)  over  all 

configurations  for  which  the  first  scatterer  is  held  fixed  at 
r^  with  a value  a-^  for  its  scattering  parameter  is  denoted  by 

^f(r,t|  q^)^>  and  defined  by 

^f(?,t|  q^)^®  J*  • • *J*  qp  • • . ,qN)p(q2» • • • »^n!  ^i^d^2*  * 'd^N» 

and  so  on. 

If  it  is  assumed  that  the  locations  and  radii  of  the 
scatterers  are  statistically  independent,  then  by  definition  we 
have 

p(q\,..-,qN)  ■ p(q\)  p(q2)---  p(^n)  * <18) 

that  is,  the  joint  probability  density  function  is  expressible 
as  the  product  of  the  individual  probability  density  functions 
for  each  scatterer. 

C.  ACOUSTIC  PULSES 

Let  Y(r, t)  be  the  acoustic  pressure  at  the  space  point 
r and  at  the  time  t.  We  will  assume  that  Y(?,t)  is  a real 
function  of  its  arguments  r and  t.  Since  we  are  interested  in 
pulse  wave  forms,  we  shall  use  the  Fourier  integral 
representation  for  'Jr(?,  t). 
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That  is,  we  let 


111 

VI 


-iojt 


*»  w jl 

(r,t)  = ^ Re  J*  i(r,iu)  e dm  , 


(19) 


where  ty(r,  tw)  is  the  Fourier  transform  of  Y(r, t)  and  is  given 

by  the  inversion  formula 

® iujt 

^(r,uu)  = J Y(r',t)  e dt  . 


(20) 


The  symbol  Re  denotes  the  real  part  of  the  integral  which  is  a 
complex  number  and  the  real  quantity  m/2TT  is  the  frequency  of 
the  sound  wave. 

The  Fourier  transform  ^(?,  m)  must  satisfy  the 
Helmholtz  equation 


y2  ♦(?,««)  + k2  *(?»«>)  * o , 


(21) 


2 11 

where y is  the  Laplacian  operator.  If  absorption  is  present, 
then  k is  a complex  number,  say 


k «*  K+ia  , 


(22) 


where  a is  the  absorption  coefficient  and  tc  is  the  wave  number 
which  is  related  to  the  wavelength  \ of  the  sound  wave  by  the 
expression 


(23) 


(24) 


The  phase  velocity  Cp  of  the  sound  wave  is  then  given  by 


= — iT“  • 

P K 

In  analogy  with  Eq.  24  a complex  velocity  c can  be 


defined  by 


Upon  using  Eqs.  22  and  24,  we  can  express  c in  the  form 


(25) 


1 + i ttCp/ uu 


(26) 


In  general  the  medium  will  be  dispersive,  so  that  both  Cp  and 
a will  be  a function  of  w. 

Now  consider  a medium  in  which  there  are  embedded 
N point  scatterers  which  are  randomly  distributed  with  respect 
to  position  and  with  respect  to  the  scattering  parameter  a^. 
Let  Y(r,  t|  ^, . . . ,^)  be  a scalar  field  (acoustic  pressure 
field)  produced  by  the  multiple  scattering  of  an  incident  wave 
Y^(r,t)  by  a configuration  of  scatterers.  Let  us  take  the 
ensemble  average  of  Eq.  19. 

{v(?,t))  YC^tl^,.  • • ,\)  P(S1f*»'3N) 


"/•••/  dV**d?N  ‘ J d 

o 

Upon  interchanging  the  order  of  integration,  we  obtain 


I 

I 


i » 


I 

i 

I 

I 


/ \ i -i^t  .s  . 

\V(vyty=  ^-Re  J e t(r,wl  3L, . . . ,^N) 

o 

* p(ci^» • • • dcf^...dq^ 
or 

<(y(r,t))=  ~~Re  J ^(r,w)^  e dm  . (27) 

Consequently  in  developing  the  theory  of  multiple  scattering, 

we  need  only  be  concerned  with  the  Fourier  components 

<^i(r,o)|>.  A more  rigorous  discussion  of  the  Fourier  transforms 

12 

of  stochastic  processes  is  given  by  Papoulis 

D.  EQUATIONS  OF  MULTIPLE  SCATTERING 

Consider  a homogeneous,  isotropic,  non-absorbing 
medium  capable  of  sustaining  wave  motion  according  to  the 
Helmholtz  equation 

y2  'Kr’,uu)  + k2  ty(r , id)  = o (28) 

where 

ko  - — , (29) 

c 
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(30) 


Bad  being  the  adiabatic  bulk  modulus  and  pQ  the  mass  density 
of  the  liquid  medium. 


Suppose  there  are  embedded  in  this  medium  N point 

scatterers  which  are  randomly  distributed  with  respect  to 

position  ? and  scattering  parameter  am  (1  <m^N) . Under  the 

influence  of  an  incident  wave  ^(r*,  uu)  and  the  scattering  from 

other  scatterers,  a scattered  wave  ^ (r,tujqm^q^, . . . ,<i  ) is 

generated  by  the  mc  scatterer.  Here  the  first  argument  ? 

specifies  the  field  point  of  evaluation,  the  second  argument  w 

is  the  angular  frequency  of  the  incident  and  scattered  radiation, 

gives  the  location  and  value  of  the  scattering  parameter  of 

the  scatterer  originating  the  radiation  and  , . . . ,q  indicate 

•**  N 

the  dependence  of  the  scattered  wave  on  the  specific 
configuration  chosen.  Employing  the  radiation  condition,  the 
scattered  wave  has  the  form  of  outgoing  radiation  and  is  a 
regular  solution  of  Eq.  28  everywhere  but  at  r = r , where  a 
singularity  is  present. 

The  properties  of  a single  scatterer  are  assumed 
known,  so  that  a rule  is  available  relating  the  scattered  wave 
and  the  exciting  field  ^(r,  w)  ) acting  on  the  mth 

scatterer  to  produce  scattering.  This  rule  defines  a linear 
scattering  operator  T (^  ) by  the  relation 


%1\>- • • >\)  ■ T($m)  'K^lqJ  qi> • • • >q^) 


(31) 


18 


I 


< • 


The  exciting  field  is  assumed  to  be  a regular  solution  of  the 
Helmholtz  equation,  Eq.  28,  in  the  neighborhood  of  r^  in  order 
that  the  single  scatterer  computation  be  applicable. 

By  the  principle  of  superposition,  the  total  field 
i|((r,«)|q^, . . . ,q  ) may  now  be  written  as 

U(r>,ui|'31,  . . . ,'?N)  = ^(r,®)  +\  f (r>lq^l  qj^,.  . . ,q^)  , (32) 


m=i 


or  using  Eq.  31, 


N 


«r,a.|t1.....tN)  - ♦*(?,»)+  £ T(qJ„)  ♦E(?>“>fqj31>.  • • .tN)  .(33) 

m=l 

th 

Now  the  exciting  field  acting  on  the  m scatterer 

th 

is  just  the  total  field  minus  the  wave  scattered  by  the  m 
scatterer,  i. e. , 

A?’ , wl  qm\ cfi , . . • «*l  q”i>  • • • T?n) - ^S(?,  ?1 » ‘ • <34> 

Using  Eqs.  31  and  33,  this  last  relation  can  be  written  as 

N 

^(r,«f§m\qi,...  ,^N)“^(r,®)+  Y,  *<3*.  ) | ?i,..., qN).  (35) 

1 

TO  =1 
1 

(m  ^m) 
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Eq.  35  represents  a set  of  N algebraic  equations  for 
the  N function  ^(r,uuj  (J  | q^,  . . . » l^m^N.  Once  these  N 

functions,  representing  the  exciting  field  on  each  scatterer, 
are  known,  they  can  be  substituted  into  Eq.  33  to  give  the  total 
field. 

These  relations  account  completely  for  the  effect  on 
each. scatterer  due  to  the  presence  of  other  scatterers.  In 
terms  of  multiple  orders  of  scattering  approach  where 
primary  (or  first  order)  scattering  is  due  to  the  incident  wave 
alone,  secondary  (or  second  order)  scattering  represents  one 
rescattering  of  the  primary  waves,  and  so  on,  all  orders  are 
included  in  Eqs  33  and  35. 


To  summarize,  the  self-consistent  field  equations  of 
multiple  scattering  are  given  by 

♦(£*»  “>1^,  • • • (£,“)+  ^ T(qm)  ^(r^lq^l^i,.  . . ,tN)  , (36) 

m=l 


. . . ,‘V-tt(2,»)+^  T^m'  * l^l»*  ' * • <37> 

m*  =1 

I 

(m  ^m) 


In  order  to  gain  some  insight  into  Eqs.  36  and  37 
let  us  consider  the  case  when  N=2.  Eq.  36  becomes 
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t(£®l'31i&2)  = ^(r'^+TCqp  ^(r^jq^q^q^) 


+T(q^)  ^(r^\\\^2)  , 


(38) 


and  Eq.  37  becomes 

♦E(r>»“l"^l!  ?1  >^2>  “ ’*l(?.“1)  + T(?2)  'tE(r^u)(q2jq17q2)  » 


♦E(£>lq2;'qi,'q2>  * + T(q’)  i»E(?>!  q^qj.qp  . (40) 

Solving  Eqs.  39  and  40  for  'r(r,  uujq^  q^,q2)  and 
^(r^ujjq^  9x^2^  » yields 


^ ^ ^ ^(r^w)  + T(q’2)  ^(r’w) 


^E(r>,uj|q1|q1,q2)  = 


1 - T(q2)  T(qp 


(41) 


,E,->  t>  ~>  t(r,w)  + T(qx)  ^(r , uj) 

* (r>uj|q2|qi’,2) 1 : t(?^7t(^27 


(42) 


The  total  field  is  obtained  by  substituting  these 
values,  Eqs.  41  and  42,  into  Eq.  38. 


'Kr*,wlq1,q2)  * ^(r>)  + T(^) 


+ T(<f2) 


^i(?,U))+T(q>2)^i(^,uJ) 
1 - T(q2)  T(qx) 


I ^(r>)  + T(q1)^i(r»  < (43) 

I 1 * T(q\)  T(q2) 
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To  show  explicitly  that  Eq . 43  includes  all  orders 
of  multiple  scattering,  we  expand  the  denominator  of  the  last 
two  terms  on  the  right  side  of  Eq.  43  in  a power  series. 
Assuming  that  T(q'1)  T(q2)  = T(q2)  T(q^)  , we  then  have  for 
both  denominators 


T(q2)  T(qp  ^ = 1 + T(q2>  T(q?1)+T(q2)T(q  1)T(q2)T(q x) 


provided  that  jT(q2)  T(q’j_)!  < 1.  Thus  Eq . 43  becomes 


'Hr>|q1,q2)  - ^(r ,w)+T(q1)  ^(r,  iu)+T(q2)  ^(r.U)) 

+T(q1)  ,T(q2)  ^(r»"  +T(q2)  [ T^) 

+T(q\)  {T(q2)[T(q1)1fi (r,0))]] 

+T(q2)  {T(q1)'T(q2)'lli(?,a))’]  +.  . . , (45) 

The  term  '|rL(r,  m)  is  the  value  of  the  incident  wave  at 
the  point  r after  having  undergone  no  scatterings;  T(q^)  <}x(r,<ju) 
is  the  value  of  the  incident  wave  at  the  point  r after  having 
been  scattered  once  by  the  scatterer  located  at  r-^;  T(q2) 'ii'L (r , w) 
is  the  value  of  the  incident  wave  at  the  point  r after  having 
• been  scattered  once  by  the  scatterer  located  at  r2;  T(q^) 

x J(q2)  l1(r,eu)j  is  the  value  of  the  incident  wave  at  the  point  r 


after  having  been  scattered  twice  - first  by  the  scatterer 
located  at  ri>  and  then  by  the  scatterer  located  at  r^;  etc. 

It  is  important  to  note  that  the  explicit  resolution 
of  the  total  field  into  multiple  orders  of  scattering,  Eq.  45, 
is  valid  only  if  the  interaction  between  the  scatterers  is 
sufficiently  weak  so  that  the  condition  | T ( ) T(q2 )|  < 1 
holds.  If  the  wave  about  each  scatterer  is  profoundly  influenced 
by  the  presence  of  the  other  scatterer,  i.e.,  if  there  is  a strong 
interaction  between  the  scatterers,  then  this  condition  may  not 
be  met  and  the  explicit  resolution  into  multiple  orders  of 
scattering  loses  its  meaning.  In  this  case  the  more  rigorous 
Eq.  43  must  be  used.  While  this  equation  accounts  for  all 
orders  of  multiple  scattering,  it  does  so  implicitly. 

If  the  scattering  is  extremely  weak,  i.e.,  if 
|T(?X)  T(cf2)|  <<1,  then  all  second  and  higher  order  scattering 
terms  can  be  neglected  in  Eq.  45  and  we  have 

♦(r»‘Mlq>l»q2>"^’(^,u,)+T(^l>  ^(^“O+TCq^)  . (46) 

This  is  called  the  Bom  approximation. 

E.  AVERAGED  FIELDS 

Because  of  the  complicated  nature  of  Eqs.  36  and  37, 
it  does  not  appear  feasible  to  attempt  to  invert  Eqs.  37  to 
obtain  explicit  expression  for  the  exciting  fields.  Instead,  the 


equations  will  be  averaged  as  they  stand.  First  let  us  average 
Eq.  36  for  the  total  field: 

<*(?,««)>“  J*.  . ./’Kr.uJlqp  . . . ,qN)p(q\,  • . . f^)dq^.  . .dq^ 


‘*X  (?,»)+ 


£ J-  • .JT(qm)  '^(r.yJlq^tq^, . . . ,qN)p(q1>.  . . ,qN)dq]  . . .dqj. 


f(r,uJ)+^  J* . . • J*T  ( cfm ) ^(r^jq^j  q\,  . . . ,qN) 


x p^m)p(qi»  * * ’ qm+i- • * j qm)  dq^ . . . dq 


-?  , -> 


= ^(?,ui)+  ^jTCq^p^dqJ.  . .Jdq^.  . . dq^dq^ . . . dq^ 


x ^(r»u,i  qm!  <Jir  • • »qN)p<?i>  • • • » 1 * %i+i * * • • >%\ 


= f (r,aJ)+}_  jT(q’  )p(qm)<^E(?,«)|q  |q  )>  dq^  . 


(47) 


Here  use  has  been  made  of  Eqs.  6,  16  and  17.  Thus  to  evaluate 
the  ensemble  average  of  the  total  field,  (ty(r,u)))>,  we  need  to 
determine  the  partial  averages  ^V^(r, iu|q  | q )^for  1 ^m^N.  The 
quantity <iiE(r, tu|  q | q is  the  ensemble  average  of  the  exciting 


j i 

field  acting  on  the  nr  scatterer  when  the  mc  scatterer  is  held 
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kU>i  ■**»»»> 


fixed.  Let  us  compute  the  first  of  these  quantities,  viz., 
<ifE(r,u)|q1|q1)>  • 


From  Eq.  37,  we  have 


N 


• >qN) 


m-2 


Consequently, 


^(r.aijq^q^y)  =J\  . .J*  ( ?>  “>  | q i f q i ’ * * * >qu)P(q2»  * ' ' ,qNl  ql)dq2‘  # ,dqN 


N 


♦i(r,w)+^  J-  • •i'T(qrn^  * >qN} 

m-2 


x p(q2» • • • »qNi ql)dq2' * ,dqN 


A slight  generalization  of  Eq.  13  yields 


(48) 


p(q2»  • • • >qN|  qi)=P (qml  q\)p(q2»  • * * ’^m-l’^m+l’*  * * ’^nI^I’^  * 


Thus  Eq.  48  can  be  written  as 
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N 

V 


<(r(r,Ui|  qL|  =v1(r,aj)+^  JT(qm)p(qm|  q^dq 

m=2 


->> 

m 


->  , -o  -> 


XJ • • •/ • • • >qjq)P(^2»  ‘ * * * ' * ,C*N  I ^1’ 

Xdq2 . . -dq^dq^. • • d<3N 


<'t'E(r\oi|q\|q1)>  =^(r, u>)+  V J*T (q^) p (q^  ( q\) 

m=2 

X^E(r,U)jqm|q1>qm)>dq>m  . (50) 

Thus  we  get  the  very  unsatisfactory  result  that  the 
average  of  the  exciting  field  with  one  scatterer  fixed  is  given 
in  terms  of  the  average  of  the  exciting  field  with  two  scatterers 
fixed.  Further  calculation  shows  that  higher  partially  averaged 
exciting  fields  will  involve  the  integral  of  the  exciting  field 
with  one  additional  scatterer  held  fixed.  This  lack  of 
completeness  is  the  basic  difficulty  encountered  in  the  implicit 
approach  to  multiple  scattering.  Lax  has  suggested  breaking  off 
this  hierarchy  of  equations  at  some  stage  by  arbitrarily  replacing 
the  exciting  field  in  an  integrand  by  the  corresponding  field  with 
one  less  scatterer  held  fixed.  The  resulting  equation  is  then 
solved,  and  each  of  the  preceding  equations  solved  by  quadrature. 
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The  system  of  exact  equations  governing  multiple 
scattering  is  now  complete.  There  are  three  known  methods  for 
obtaining  a solution  to  these  equations.  First,  one  can  in 
principle  obtain  an  exact  solution  by  solving  the  implicit 
initial  representation,  Eq.  37,  for  the  N quantities 
^(r,  ou|  qm|  q^, . . . ,q^) . This  method  was  illustrated  by  the 
example  for'  the  case  when  N = 2 with  the  solutions  being  given 
by  Eqs.  41  and  42.  The  solutions  for  ^(r,  rnj  q^|  q3^,  . . .,q  ) , 
obtained  by  solving  the  system  of  equations,  Eq.  37, » can  then  be 
substituted  into  Eq.  47  to  give  an  exact  expression  for  the 
ensemble  average  of  the  total  field.  However,  this  scheme  must 
be  discarded  because  the  labor  is  prohibitive.  Second,  one 
might  iterate  the  system  (eq.  37),  replacing  the  exciting  field 
in  the  summation  of  one  equation  by  the  right-hand  side  of  the 
following  equation,  and  continuing  to  replace  the  exciting  field 
wherever  it  occurs,  ultimately  obtaining  an  infinite  series 
representation  for  each  of  the  exciting  fields  which  involves 
only  operations  on  ^ ( r*,  uj)  . This  result  expresses  the  exciting 
field  in  terms  of  multiple  orders  of  scattering  and  leads  to  the 
system  of  equations 
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N 


♦E  ( f , U.  | ? | f qN  ) - i)1  (?,»)+ 1_  T 1 ' ) (f1  ( ?,  ® ) 

m'=l 

(mVm) 


N 


N 


+Z  I(C>Z  T(5^' )«*■(?, »)+...  . 

m"-l  m'=l 

(m'Vm)  (mVni) 


(51) 


This  result  is  then  substituted  into  Eq.  47  to  obtain  the  average 
for  the  total  field.  Such  series  rapidly  become  cumbersome  to 
treat  exactly  when  more  than  one  or  two  terms  are  considered, 
so  that  this  scheme  is  only  feasible  when  multiple  scattering 
effects  are  very  weak.  Finally,  one  can  average  Eq . 37  as  it 
stands.  This  is  the  technique  that  was  used  in  this  section  and 
it  was  seen  to  lead  to  a complicated  hierarchy  of  equations. 

Even  if  Lax's  suggestion  is  used,  the  labor  may  be  prohibitive. 
Thus,  if  we  are  to  arrive  at  a solution,  some  kind  of  approximation 
must  be  used.  One  such  approximation  will  be  discussed  in  the 
next  section. 


F.  APPROXIMATE  EQUATIONS 

We  begin  by  expressing  the  exciting  field 
irX?, <ju]q^|  q-^, . . . ,q^)incident  on  the  mc  scatterer  in  terms  of  the 
total  field  ^'(r,  uu | , . . . » ]_ » ?m+]_ » • • • »?N)  that  would  exist  in  the 
neighborhood  of  the  mc  scatterer  if  the  mL  scatterer  were  not 
present.  That  is,  we  first  remove  the  in  scatterer  from  our 
configuration.  Then  let  t(r»wJq\, . . . • • • ><?N)  he  the 
total  field  in  the  neighborhood  of  the  position  where  the 
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th 


m scatterer  was  originally  before  removal.  We  now  insert  the 


th 


m scatterer  back  into  its  original  position.  Due  to  its 


presence,  the  wave  '['(r,  uu|  q^ , . . . • • • >9^)  * which  existed 


th 


in  the  vacant  neighborhood  of  the  m scatterer,  is  scattered 


th 


by  the  m scatterer  and  then  rescattered  by  the  other 


scatterers,  consquently,  changing  the  value  of  the  wave 

th  - . . -<>  ->  -o 


incident  on  the  in  scatterer  from  ^(r,  w|q^, . . . ,qm_^  qm+^,  * * *,9^) 
to  ^(r,  w|q  | q^, . . . ,q^)  • We  can  express  this  as 


^E(r,  U)|qml  qp  . . . ,qN)  = 'Kr,<u|q1,  . . . VkL>  . . . ,qN> 


N 


+ > 


XS(r,a)lqk|q’1,...,qN)  , 


(52) 


k=l 

(k/m) 


where  X (r, uu|  qk|  q^,  . . . ,qN)  is  the  wave  scattered  by  the  k^ 
scatterer.  Let  x*  (r,  w|  cfk)  q^, . . . ,qN)  be  the  wave  incident  on 


.th 


the  k scatterer  to  produce  the  scattered  wave 


X (£">  ^1  9k|  qp  • • • >9jq)  > i.e.,  let 


-■> 


X (r,w|qk|q1,...,qN)=T(qk)X  (r,  w|  qR|  qp . . . ,qN> 


(53) 


From  the  above  discussion  we  can  write 

E / r?  ...i?  i -•>  n m,r>  \ -*>  r?  -■> 


X (r>ai9k|qp«»*  >?jq)~T(qm)  *(r,w|q1,  . . >%H-1 » • • • »9jj) 


N 


+ } T(q\  ) xE( r , w|  q j | q\ , . . . , qN)  . 

(j#n,k) 


(54) 


'i 


1 


1 


j 
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We  now  apply  the  same  iteration  technique  to  Eq . 54  that  we 
used  to  obtain  Eq.  51  from  Eq . 37.  This  yields 


XE(r,U)|q>kJq1,  • • • ,?N)“T(qm)  V-l’^m+l’  * 


N 

+L  T(?j>T(0  ^?»a)1^i %i-i’%»+i’ 

j=l 


+ 


(55) 


If  we  now  combine  Eqs.  52,  53  and  55,  we  have 
, u)j  qmS  qf,  . . . ,qN>“  K?»  • • • »?m-l»^m+l*  * * ' 


+ V T(q.  ) 1 + T(q' ) + ■ T(?,)  ..  . 

k^m  k u j^k,m  J J5^k,m  J n^j 


j^k,m  J ,m 

+ ^ T<&)  1 T(q’)  ^ T(q’ )+...] 

j^k,m  J n^j,m  n p^n,m 


T(?n) 


X T(qm)  |q\,  • • • »?m-l>?m+l»  * • * >?N) 


(56) 


It  should  be  stressed  that  this  representation,  while  not  in 
closed  form,  is  exact. 
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Now  a great  simplification  in  the  equations  of 
multiple  scattering  would  result  if  we  could  make  the 
approximation 

tu|  qm|  q^, . . . ,qN>  * K?,  «>1  » • • • Vl-1*  * • ‘ »^N>  ’ (57> 

obtained  by  neglecting  all  but  the  first  term  on  the  right  side 
of  Eq.  56.  Unfortunately,  there  are  no  known  general  conditions 
which  tell  us  when  these  terms  can  be  neglected.  However,  for 
the  case  of  isotropic  scattering  an  estimate  can  be  made  and  this 
will  be  discussed  in  Section  G. 

Using  the  approximation  given  by  Eq.  57,  that  the 
exciting  field  on  a scatterer  may  be  approximated  by  the  total 
field  existent  in  the  neighborhood  if  the  scatterer  were  not 
there,  we  can  write  Eq.  36  as 


K?>lq\,  • • • ,qfl)  =^X(r ,id)+V  T(qm) 

m=l 


x ^(r,  w| q^, . . . J^m+l » • • • 


-■5 


(58) 


Now  let  us  assume  that  the  locations  and  radii  of 
the  scatterers  are  statistically  independent  so  that  Eq.  18 
holds.  Then  taking  the  ensemble  average  of  Eq.  58  yields 

0(?,,x’)>=  **<?.«>)  + /P(%i>T(^m>d%i  J • • -dVldi5i+r  * * d^N 

m=l 

x P(q^)  • • • P(<5m-i)P(^m+l^  ’ ’ *P^N^  ^i»  • • • ’^m-l’^m+l’  * * * ’^N^  * 
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The  (N-l)  - fold  integral  in  Eq.  59  is  independent  of  and, 
in  fact,  is  precisely  the  average  total  field  for  the  same 
problem  with  one  less  scatterer  involved.  As  the  number  of 
scatterers  increases,  this  quantity  must  differ  from  the  total 
field  with  N scatterers  by  terms  of  order  1/N.  Thus,  assuming 
a large  number  of  scatterers,  we  make  the  approximation 

^'^(r , uj)}  = J . . . J* dq^ . . . dcf^  s'(r,u>|q^,  . . . , q^)p (q^)  • • • P (^Ijj) 

=J*  • -Jd9]_-  • • d?m_ ld%i+l * * •d^N^r,UJ^l’  ‘ - ,qm-l’qnH-i’  - * * ,qN^ 


l. 
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x p(q?1)...p(q^1_1)p(qJ1+1)...p(qN)  . (60) 


Consequently,  Eq.  59  becomes 
N 

w)'/  = ^(r,u)>  + ) jT(qm)p(q>m)<K?,i«)>dqr] 


m 


m=l 


Now  using  Eq.  5 and  assuming  that 


p(qi)  = p(q'2)  = •••=  » 

n(q^)  = n(q2 ) ~ • ■ • = ^(q^q)  » 


We  can  write  Eq.  61  as 

<t(r,« >)>-  ^(r,^)  + jT(q  )n(f  )<*(r\u>))dtf  . 
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(61) 


(62) 


(63) 


This  integral  equation  will  be  the  starting  point  for  our 
investigation  of  the  multiple  scattering  from  fish  schools, 


G.  ISOTROPIC  SCATTERING 


For  the  case  of  isotropic  scattering,  the  operator 


T(q  ) is  defined  by 
m 


t(%>  (E(?,  »l  %i  q\ ?N)=s  (4, »)  “ 1^,1  fi q»)f= 


itc  1 r-r  I 
o 1 m1 


. (64) 


Thus,  the  integrand  on  the  right  side  of  Eq.  63  becomes 


t(!  )(Hr, oj))=t(^  )J . . ./♦(?» «>|qi» • • • »qN)p(q\, • • . .qpdq^. . .dq>N 


"J • . • • • >qN)T(q  ) K?,uj|qp  . . . »qN)dqr  . .dq^ 


i K i r-  r I 


-/••■J,p(q1»...»qN)g(q  ,(«|q1,  • • . ,qN)  e ° dq^ . . . dqN 


i i*  !r-r  | i 

= g(q  e (K*  ,w)>  • 

I r-r  I 


The  operator  T(c?  ) commutes  with  the  N-fold  integration  since 
cf  is  independent  of  q^,...,q^  and,  so  far  as  the  N-fold 
integration  is  concerned,  the  variables  q^,...,q^  are  simply 
dummy  variables. 


(65) 
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Consequently,  using  Eq . 65,  our  integral  equation 
for  multiple  scattering,  Eq.  63,  becomes  for  the  case  of 
isotropic  scattering 


| ? ? | 

<Kr»>=  '^(r\«0+J7n(?  >a  )S(?  >a  dr'da'.  (66) 

Jr-r  | 


Defining 


1 

•7  . \ / ' 


G(?  ,iu)=Jn(r  ,a  )g(r  ,a  ,x)da  , 


(67) 


Eq.  66  takes  the  final  form 


r , uj)>  = r ( ?,  a ) +J G ( r » <(  K r ',  «))  e„_!!.r_ 

|r-'r| 


iK  _|  r-r  I i 


dr  . 


To  see  the  physical  significance  of  this  equation. 


let  us  apply  the  operator 


(68) 


2 . ..2 


V + Kc 


(69) 


to  both  sides  of  the  integral  equation,  Eq.  68,  remembering 


that 


2 2 lK0lr"r  I 

<V  + V>  *».(M  ) , 


r-r 


where  6(r-?  ) is  the  three  dimensional  Dirac  6-function 
defined  by  the  equation 


(70) 
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J*  f(r  )6(r-r  )dr  = f(r)  , if  rev 
V 

= o,  if  r i V . (71) 

We  then  obtain 

y2<((r,m)>  + k2(r,u>)0(r,w)>=  o , (72) 

where 

k2(r,m)  = k2  +4rt  G(r, uj)  . (73) 

In  general,  g(r,a,u)  will  be  complex,  so  that  k(r,w)  will  be 
complex.  Thus,  wc  see  that  satisfies  the  Helmholtz 

equation  having  a complex  wave  number  that  depends,  in  general, 
on  position  and  frequency.  Consequently,  the  incident  wave  and 
the  scattered  waves  from  all  the  scatterers  interfere,  on  the 
average,  to  form  a new  wave  traveling  at  a different  phase 
velocity  and  undergoing  attenuation.  This  wave  will  display  the 
reflection  and  refraction  aspects  of  coherent  scattering  at 
surfaces  of  discontinuity. 

The  problem  of  finding  the  average  value  of  the  wave 
function  has  been  essentially  reduced  to  solving  a boundary 
value  problem  in  the  wave  equation.  The  boundary  conditions  are 
implied  in  the  integral  equation  itself  and  depend  on  the 
function  G(?,  uu)  . If  G(r, u)  is  everywhere  continuous  and 
approaches  a constant  value  or  zero  at  infinity,  then  the 


1 


: 


\ 

I 


boundary  conditions  are  that  ^(r,  iu)')>-  be  everywhere 

continuous  and  have  a continuous  gradient  and  at  infinity 
represent  outward  traveling  waves.  in  another  important 
case,  that  in  which  G(r,  uo)  is  sectionally  continuous,  the 
boundary  conditions  are  that  ^(r , uu))>  - ^(r,  iu)  be  everywhere 
continuous  and  have  a continuous  normal  derivative  across  a 
surface  of  discontinuity  of  G(r,  id)  and  at  infinity  represent 
outward  traveling  waves,  provided  again  that  G(r,  id)  approaches 
a constant  value  or  zero  at  infinity.  In  both  cases,  of  course, 
<♦(?  , «j)^>  must  approach  ^(r,  ud)  as  G(r,w)  becomes  zero  everywhere. 

We  note  also  that  the  integral  equation,  Eq.  68,  can 
in  principle  be  solved  directly.  One  method  of  solution  is  the 
Liouville-Neumann  method  of  successive  substitutions  or  iteration 
method.  It  consists  in  repeatedly  substituting  the  expression 
for  u )')  as  given  by  the  right  side  of  the  integral  equation 
for  ^ty(r  >w))>  under  the  integral  sign,  thus  yielding  in  our 
case  the  infinite  series 

CO 

<Kr,  w)>=  ^ ^ (?><*>)  » (74) 

m=o 


where 


JG(  r,a>) 


m 


-i(r> 


m) 


\r-r  | 


iK  | r-r  | 

n * 


dr  , (m^o) 


(75) 


which,  if  it  converges  uniformly,  is  the  desired  solution. 
In  the  above  notation  (r,  w)  = w) . If  the  scattering 

is  very  weak,  it  may  be  possible  to  neglect  all  terms  for 
which  m > 2.  In  this  case  Eq . 74  takes  the  form 


<^(r,uu)>  « ^(r,ai)  + jG(r,ai)  ^(r,^) 


i K | r - r | 
e ° 

— T-Ti dr  . 

I ?-r  | 


(76) 


This  approximation  is  known  as  the  Bom  approximation. 

In  order  to  obtain  the  approximate  integral 
equation,  Eq.  63,  or  for  the  case  of  isotropic  scattering,  Eq.  68, 
we  made  the  assumption  that 

^(?,wl$nl*L,---»HN>  ~ n?,w|qls . • • ,qm-i,qm+i> • • • ,qN)  . 

Let  us  now  examine  the  conditions  under  which  this 
approximation  is  valid  for  the  case  of  isotropic  scattering. 

The  argument  while  plausible  is  by  no  means  rigorous.  Consider 
the  first  term,  say  cp  , of  the  series  in  Eq.  56,  i.e.,  let 


cpG(?)  = T(?k>  <77> 

k^m 

Each  term  in  the  summation  represents  rescattering  from  one  of 

the  obstacles  of  an  outgoing  scattered  wave  from  r^.  Furthermore, 

since  the  exciting  field  is  only  required  in  the  immediate 

neighborhood  of  the  scatterer,  in  question,  in  this  case  the 
til 

m according  to  Eq.  56,  we  need  only  estimate  cpQ(?)  in  this 
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neighborhood;  for  simplicity  consider  cp0(r^).  Surfaces  of  equal 
phase  of  the  outgoing  wave  will  be  spherical  so  that  the  phase 
of  rescattered  waves  returning  to  r retard  according  to  the  round 
trip  distance  2 [ r ^ from  to  the  individual  scatterer.  A 

set  of  concentric  half-period  zones  may  be  constructed  about 
r^,  each  zone  defined  by  the  requirement  that  its  rescattered 
waves  are  no  more  than  half  a period  out  of  phase  at  r^.  These 
zones  are  illustrated  in  Fig.  1. 
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Suppose  there  is  a uniform  random  distribution  of 


density  n and  further  that  the  scattering  amplitude 
g(r^,a^,m),  defined  by  Eq.  64,  is  independent  of  the  position 
r^  of  the  ku  scatterer  and  that  g(a,w)  is  identical  for  all 
scatterers.  Then  using  the  fact  that  the  scattering  operator  T 
is  defined  by  Eq . 64  for  the  case  of  isotropic  scattering  and 

using  the  abbreviation  = ' * * ’%i-l’%W-l 9^)’ 

the  expression  for  cs  (r  ) becomes 


«P. 


(r  ) = 

i v m7  — 

k/m 


2iK 


g2  ><?  ) 


r,  -r 
k m 


r,  -r 
k m 


(78) 


The  contributions  from  succeeding  zones  are  the 
same  order  of  magnitude  but  alternating  in  sign,  the  square 
law  increase  of  number  of  scatterers  with  distance  being 
precisely  compensated  by  the  inverse  behavior  of  both  scattered 
and  rescattered  waves.  One  concludes  that  summation  over  the 
first  zone  should  suffice  to  give  an  order  of  magnitude  estimate 
of  cPqC?;),  and  further  approximating  the  sum  by  an  integral  we 
obtain 
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Thus  cpo(?m)  may  be  neglected  in  comparison  with 
‘ ‘ * >qm-l,qm+l’  * ’ * >qV  in  Eq*  56  Provided 
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Since  4rr|g|  can  be  shown  to  be  just  the  scattering 
cross  section  o of  a single  scatterer,  the  above  criterion 
simplifies  to 

4r~  « 1.  (81) 


The  above  estimates  have  only  been  concerned  with  the 
primary  order  of  the  multiple  orders  of  scattering  correction 
to  the  exciting  field,  i.e.,  the  first  term  of  the  infinite 


40 


it 


series  appearing  on  the  right  side  of  Eq.  56.  One  can  then 
infer  that  the  estimate  remains  valid  when  the  entire  series 
is  considered^. 


Part  II.  SCATTERING  FROM  FISH  SCHOOLS 


A.  INTRODUCTION 

The  general  theory  of  the  multiple  scattering  of 
a pulse  by  a random  collection  of  point  scatterers  that  was 
developed  in  Part  I is  now  applied  to  the  case  where  the 
scatterer  is  a bladder  fish.  The  reason  for  confining  our 
attention  to  bladder  fish  will  be  explained  later  in  this 
section.  In  general,  we  may  distinguish  between  "scattering 
groups"  and  "scattering  layers".  By  a "scattering  group"  we 
shall  mean  a collection  of  fish  such  that  the  horizontal 
dimensions  are  at  most  only  a few  times  larger  than  the  vertical 
dimension.  By  a "scattering  layer"  we  shall  mean  a collection 
of  fish  such  that  the  horizontal  dimensions  are  a great  many 
times  larger  than  the  vertical  dimension.  A "scattering  layer" 
can  be  represented  mathematically  by  a volume  bounded  by  two 
horizontal  planes  which  extend  to  infinity.  In  general,  the 
planes  will  have  an  uneven  surface.  As  far  as  this  report  is 
concerned,  a school  of  fish  will  mean  a "scattering  group" 
while  the  only  "scattering  layer"  that  we  shall  be  interested 
in  is  the  Deep  Scattering  Layer  (DSL).  Even  though  we 
mentioned  in  the  Introduction  to  this  report  that  the  DSL  does 
not  appear  to  be  a likely  source  of  false  contact  for  operation 
in  the  BB  mode  of  propagation  since  it  would  be  detected 
simultaneously  on  more  than  one  beam,  we  are  considering  it 
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briefly  in  this  report  because  of  its  importance  in  producing 
reverberation  and  because  it  is  an  immediate  consequence  of  our 
theory. 

Large  areas  of  the  deep  ocean  have  a virtually 

i 

continuous  DSL.  These  are  formed  of  various  sorts  of  marine 
organism,  but  especially  of  bathypelagic  fish.  Out  of  the  many 
studies  attention  should  be  drawn  to  those  that  have  looked  at 
the  frequency  structure,  and  by  using  wide-band  underwater 
explosion  sources  have  demonstrated  resonances.  * * These 

results  are  in  fact  the  strongest  evidence  that  resonant  swim 
bladders  do  play  an  important  part  in  fish  acoustics.  Extensive 
scattering  layers  can  also  occur  in  shallow  coastal  waters  and 
are  frequently  associated  with  summer  thermoclines . 

Tightly  packed  schools  of  fish  can  occur  in  both  the 
deep  ocean  and  in  shallow  water,  and  they  come  in  a great 

variety  of  shapes,  sizes  and  concentrations.  The  records  , .> 

obtained  by  Voglis  and  Cook^  show  some  of  this  variety  for 

' 

schools  of  pilchard.  The  school  sizes  which  have  been  seen  with 
the  Voglis  equipment  vary  from  a few  yards  across  to  those  with 
a maximum  dimension  of  about  40  yards  (the  soundings  were  taken 
looking  down  obliquely  from  above).  The  shapes  vary  from  the 
near  circular  to  the  line,  a surprising  number  have  U-shape. 

The  differences  in  the  apparent  strength  of  return,  as  measured 
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by  Voglis,  must  not  be  interpreted  solely  in  terms  of  mean  fish 
concentration,  it  may  also  depend  on  the  roughness  or 
diffuseness  of  the  school  boundary. 

While  scattering  of  sound  by  fish  is  often  observed, 
the  possibility  of  significant  attenuation  due  to  fish  has  not 
been  discussed  in  the  literature.  The  status  of  the  subject  of 
low  frequency  attenuation  is  that  in  deep  water  its  high  value 
is  still  unexplained,  and  in  shallow  water  Weston  and  his 
colleaques  have  made  direct  observations  proving  that  a large 
part  of  the  long  range  attenuation  is  sometimes  due  to  fish. 

While  this  report  is  only  concerned  with  calculating  the 
first  order  statistics  of  the  signal  scattered  by  a fish  school, 
future  work  will  deal  with  the  second  order  statistics  of  the 
scattered  signal.  Consequently,  once  the  second  order  statistics 
are  known,  the  question  of  whether  there  is  significant 
attentuation  due  to  fish  schools  will  be  answered. 

Now  let  us  examine  the  theoretical  and  experimental 
work  which  has  been  done  on  sound  scattering  from  fish.  For 
the  deep  scattering  layer,  we  have  already  mentioned  that 
sound  scattering  measurements  have  been  made,  e.g.,  by 
Hersey,  et.  al.}^  Marshall,  et.  al.}^  and  Andreeva^  The 
scattering  of  a pulse  from  the  DSL,  neglecting  multiple 
scattering,  has  been  theoretically  calculated  by  Glotov^ 
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Glotov,  et.  al.^  and  Kuryanov  There  seems  to  have 

been  no  work  done  which  describes  the  multiple  scattering  of 
either  an  incident  monochromatic  wave  or  an  incident  pulse  by 
the  DSL. 

No  theoretical  work  of  any  kind  appears  to  have  been 

done  on  sound  scattering  from  fish  schools  and  only  a very 

limited  amount  of  experimental  work.  Only  three  papers  could 

20 

be  found  which  reported  measured  values:  Wickham,  et.  al.  , 

21  22 

Voglis,  et.  al.  , and  McCartney,  et.  al.  Of  these  three 

papers,  only  the  first  two  and  another  paper  by  McFarland, 

23 

et.  al.  , contained  any  useful  information  on  the  structure  of 
fish  schools  and  even  this  information  was  rather  sketchy. 

The  work  which  has  been  done  on  the  scattering  of 
sound  by  a single  fish  is  considerably  more  extensive.  Most  of 
the  measurements  deal  with  the  scattering  cross-section  of 
the  fish  which  gives  information  as  to  the  form  of  the  scattering 
amplitude  g(?,a,  undefined  by  Eq.  64.  However,  there  is 
considerable  discrepancy  and  variability  among  the  results 
reported  by  the  various  investigators  which  warrants  a close 
look  at  these  papers. 

2 

First,  consider  the  work  of  Cushing,  et  al  . In  this 
series  of  experiments  the  scattering  cross  section  of  fish 
were  measured.  The  measurements  were  made  at  30  kHz  on  fish 
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which  were  between  15  cm  and  100  cm  in  length.  This  data  allows 
one  to  plot  a single  curve  of  the  scattering  cross-section  versus 
the  length  of  the  fish,  L,  for  15^L<  100  cm.  However,  Cushing 
first  extrapolates  this  curve  to  include  fish  in  the  range 
0<L<100,000  cm.  Then  he  extrapolates  in  the  frequency  domain 
by  plotting  curves  of  the  scattering  cross-section  versus  fish 
length  for  0<L<  100,000  cm  and  for  a frequency  range  from 
1 kHz  to  10  MHz.  All  this  from  one  series  of  measurements  of 
30  kHz  on  fish  between  15  cm  and  100  cm  in  length!  Moreover, 
there  is  considerable  variability  in  the  scattering  cross-section 
of  fish  of  similar  size.  For  example,  the  cross-section,  o, 
from  some  measurements  on  perch  is  given  by: 


L(cm) 

a (cm 

24.0 

3.1 

22.5 

28.7 

21.5 

71.0 

20.5 

9.8 

To  quote  Cushing:  "No  satisfactory  explanation  could  be  found 
to  account  for  these  results". 

The  experimental  work  of  Has  let  t^  is  much  more 
carefully  performed  than  that  of  Cushing.  Haslett  made 
measurements  on  guppies  and  sticklebacks  (0.935L<6.2  cm)  at 
frequencies  of  360  kHz.,  625  kHz  and  1.48  MHz.  However,  like 


Cushing,  he  also  extrapolated  his  measured  values. 

An  examination  of  the  work  of  Cushing  and  Haslett 
shows  that  the  actual  measured  values  were  made  for  fish 
lengths  L>Xo,  where  XQ  is  the  wavelength  in  water,  and  in  both 
papers  the  results  were  extrapolated  into  the  region  L<Xo. 
However,  measurements  have  been  made  in  the  region  L<Xq  by 
Coate  J and  Andreeva  and  they  disagree  with  the  extrapolated 
values  of  Cushing  and  Haslett.  Thus  it  must  be  concluded  that 
the  extrapolated  curves  of  Cushing  and  Haslett  do  not  hold 
in  the  region  L<XQ. 

Now  two  questions  must  be  asked:  (1)  What  causes  the 
difference  in  the  scattering  properties  of  fish  in  these  two 
regions;  (2)  Is  L^Xq  a valid  criterion  for  the  transition  point 
between  these  two  regions?  Unfortunately,  the  answers  to  these 
two  questions  are  not  forthcoming  and  one  can  only  speculate. 
Weston's  opinion  is  that  the  return  from  the  fish  tissue 
becomes  important  when  the  fish  dimensions  are  comparable  to 
the  wavelength,  but  in  fact  the  height  and  width  dimensions 
are  also  significant  and  the  transition  really  occupies  a wide 
frequency  interval.  Above  this  transition,  wave  interference 
effects  occur  between  different  parts  of  the  fish,  different 
parts  of  the  bladder,  and  between  fish  and  bladder.  Above 
the  transition  the  tissue  and  the  bladder  returns  have  the  same 


order  of  magnitude.  Below  the  transition  the  bladder  echo 
predominates  and  below  bladder  resonance  both  bladder  and  tissue 
are  in  the  Rayleigh  scattering  region  (cr  oc  1/X^).  Weston 
concludes  by  stating  that  "at  the  low  frequencies  (i.e.,  LAq<1) 
there  is  a great  need  for  comprehensive  target  strength 
measurements,  especially  those  covering  a wide  frequency  range 
about  the  resonance."  It  would  be  best  to  use  undamaged  live 
fish,  removing  uncertainties  such  as  the  possibility  of  bladder 
shrinkage."  This  latter  remark  was  made  since  measurements  to 
date  have  been  performed  on  dead  fish:  some  freshly  killed, 
some  kept  in  formaldehyde  for  several  weeks  before  use,  some 
kept  frozen  for  several  weeks  and  thawed  before  use  and  others 
fitted  with  artificial  swim  bladders. 

Weston's  remarks  are  borne  out  by  the  work  of  Coate 
in  the  region  L<\q.  Coate  showed  conclusively  that  for  fish 
of  length  20.5  cm  and  21.7  cm  and  for  the  frequency  range  of. 

150  Hz  to  1000  Hz,  the  scattering  was  due  entirely  to  the  air 
bladder  and  moreover,  for  purposes  of  calculation,  the  actual 
elongated  bladder  could  be  replaced  by  an  equivalent  spherical 
bladder. 

The  effect  that  these  results  have  on  the 
investigation  of  the  long  range  classification  problem  posed 
by  fish  schools  can  be  summarized  as  follows:  The  theory  of 


multiple  scattering  is  based  ultimately  on  the  scattering 
properties  of  a single  fish.  Thus  it  is  imperative  to  have  a 
reliable  model  for  single  fish  scattering.  For  definiteness 
consider  a sonar  operating  at  a frequency  of  3000  Hz  so  that 
X = 50  cm.  Using  the  approximate  criterion  above,  it  can  be 
said  that  the  scattering  from  fish  of  length  L < 50  cm  is  due 
entirely  to  the  bladder.  A detailed  mathematical  model 
describing  scattering  by  a bladder  has  been  developed  and 
agrees  quite  well  with  experimental  results.  For  fish  of 
length  L > 50  cm  the  scattering  from  the  fish  tissue  is  as 
important  as  that  from  the  bladder.  However,  no  mathematical 
model  has  been  developed  to  describe  the  scattering  process  in 
this  region.  This  is  unfortunate  since  these  larger  fish  may 
cause  more  of  a false  alarm  problem  than  the  smaller  fish. 

The  reason  for  this  is  that  experimental  results  indicate  that 
the  scattering  cross-section  increases  as  the  length  of  the 
fish  increases.  However,  it  is  not  known  if  these  larger  fish 
school  and  if  they  do,  just  how  large  the  schools  actually  are. 

Thus,  for  the  reasons  just  explained,  we  will 
consider  only  the  case  for  scattering  from  bladder  fish  in  the 
region  L/Xq  < 1.  Two  geometries  will  be  considered  for  the 
boundary  surface  of  the  fish  school:  (1)  the  fish  are  contained 
within  a spherical  region;  (2)  the  fish  are  contained  within  a 
layer  with  infinite  plane  boundaries.  In  both  cases  it  will  be 
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assumed  that  the  boundary  surfaces  are  smooth.  These  two 
geometries  are  considered  first  simply  because  they  are  the 
most  amenable  to  solution  and  understanding  the  scattering 
characteristics  of  obstacles  with  simple  shapes  and  simple 
acoustic  properties,  such  as  the  sphere  and  the  plane,  is  a 
necessary  stepping  stone  to  an  adequate  understanding  of 
scattering  from  obstacles  more  complicated  in  shape,  composition 
and  structure. 

B.  SCATTERING  FROM  A SINGLE  FISH 

In  the  region  lAq  <1,  scattering  is  due  entirely 

to  the  air  bladder  of  the  fish  and,  to  a good  approximation, 

the  actual  elongated  bladder  can  be  replaced  by  an  equivalent 

25 

spherical  air  bladder  for  purposes  of  computation  . The 
scattering  characteristics  of  the  bladder  are  completely 
specified  once  the  scattering  amplitude  g(r,a,u)),  defined  by 
Eq.  64,  is  known.  This  function  is  given  by  Weston  and 
Andreeva^  as 


g(a,iu)  = 


- 1 -i 


(82) 


where  a is  the  radius  of  the  equivalent  spherical  bladder, 
uj/2tt  is  the  frequency  of  the  incident  sound  wave,  uj^./ 2tt  is  the 
resonant  frequency  of  the  bladder  and  Q is  a damping  factor 
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which  describes  the  sharpness  of  resonance.  Note  that  we  are 
assuming  that  the  scattering  amplitude  g(a,iu)  is  independent 
of  the  position  r* of  the  fish.  Moreover,  we  will  assume  that 
the  scattering  amplitude  is  identical  for  each  fish. 

The  resonant  frequency  u)^  of  the  swim  bladder  is  given 


a,  -I. 
r a 


3yPh+4u- 


(83) 


where  PQ  is  the  density  of  water,  is  the  hydrostatic  pressure 
on  the  bladder  and  y = Cp/Cv,  is  the  ratio  of  the  specific  heat 
at  constant  pressure,  C , to  the  specific  heat  at  constant 

r 

volume,  Cv,  of  the  gas  in  the  bladder.  The  complex  shear 
modulus,  )j,  of  the  bodily  tissues  of  the  fish  which  surround 
the  air  bladder  is  written  in  terms  of  its  real  and  imaginary 
parts  as  |i  = ^(l+i^)* 

The  Q-factor,  which  describes  the  energy  losses 
incurred  in  the  pulsations  of  the  bladder,  can  be  expressed 
as  a sum  of  three  terms 


(84) 


where  Qr^  describes  energy  losses  due  to  reradiation, 
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describes  thermal  losses  and  describes  viscous  losses  in 
the  fish  tissue. 

For  radiation  damping 


a 


(85) 


For  thermal  damping 


9Ph(V-l) 

iij-g 

p m a 
Ko  r 


YV 


2m_ 


(86) 


by 


For  viscous  damping 

1 _ 4m1u2 
Q7  n "22  ' 

f po  ra 


(87) 


Here  cq  is  the  speed  of  sound  in  water  and  is  given 


c =, 


B 


ad 


'o  V po  ’ 


(88) 


where  Ba(j  is  the  adiabatic  bulk  modulus  of  water,  and  v1  is  the 
thermal  diffusivity  of  the  gas  in  the  bladder  and  is  defined  by 
, K 

v = — , (89) 


P C 
8 v 


where  Kfc  is  the  thermal  conductivity  of  the  gas  and  is  the 
density  of  the  gas  in  the  bladder. 


T 


I 

T 

C.  SCATTERING  FROM  A SPHERICAL  SCHOOL  OF  FISH 

Consider  a random  collection  of  fish,  all  with 
identical  scattering  properties,  contained  within  a spherical 
volume  of  radius  A.  Let  the  center  of  the  sphere  be  located 
at  the  origin  of  the  polar  coordinate  system  r,0,cp  (cf.  Fig.  2). 

1 • 

Since  the  fish  are  represented  acoustically  by  their  air  bladders, 
the  density  p£  of  the  medium  enclosed  by  the  spherical  boundary 
will  be  simply  the  density  of  the  mixture  of  the  water  and  the 
air  within  the  bladders.  Further,  we  assume  the  water  to  be 


nonabsorbing  and  infinite  in  extent.  Let  pQ  denote  the 
constant  density  of  the  water  and  cQ  the  constant  speed  of 
sound  in  the  water  as  defined  by  Eq.  88. 

A plane  acoustic  CW-pulse  of  frequency  w.  /2TT,  pulse 
length  2t  and  acoustic  pressure  'i'1(?,t),  traveling  parallel  to 
the  polar  axis  in  the  positive  z-direction,  impinges  upon  the 
spherical  school  of  fish.  (Choosing  the  incident  wave  in  such 
a manner  eliminates  the  dependence  on  tp . ) The  incident  pulse 
gives  rise  to  an  average  internal  wave  /y^(r,  t)^>  and  an  average 
external  scattered  wave  ^Ys(r,t)^>. 

The  incident  CW-pulse  can  be  expressed  in  the  form 


t * z/c  ' 

t (r, t)  = n cos  (Ktz  - uutt)  , 


(90) 


where  we  have  put 


Kt  " 


(91) 


The  function  Il(§)  is  defined  by 

n (?)  - i , if  |?l  < 1/2  , 

- 0 , if  |?|  > 1/2. 

Now  let  us  compute  the  Fourier  transform  of 
¥*"(?, t).  Using  Eq.  20,  we  obtain 


(92) 
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*• 


00  . iu>t 

^(r,  u')  = J ^(r,  t)  e dt 

— CD 

t - z/c 


2 r 


lUJt 

COS  (KtZ  - u’tz)  e dt 


Z/C0+T 


(r,ui)  = J cos  (K^z  - 


iojt 

uutz)  e dt 


z/cq-t 


. . i K Z 

'^(r,  a1)  = e 0 . 


(96) 


Let  n(x?'  ,a')  da'  be  the  average  number  of  scatterers 
per  unit  volume  in  the  neighborhood  of  the  point  r'  having 
their  scattering  parameters  lying  between  a'  and  a'+da'.  The 
function  n(r?,,a')  is  just  the  number  density  function 
introduced  in  Part  I,  Section  B.  Then  using  Eq.  82  for  the 
scattering  amplitude  of  a single  fish  g(a',  to  ),  the  scattering 
amplitude  for  the  school  G(r',u'),  defined  by  Eq.  67,  becomes 


® n(r 1 ,a  * ) a 1 da 1 

— : — r 

° K-J-WO 

V / i “Q  / 


(97) 


We  will  assume  that  there  are  no  fish  outside  the 
spherical  volume  of  radius  A and  that  the  fish  inside  the 
spherical  volume  are  uniformly  distributed,  i.e.,  we  put 


n(^,,a')  «n6(a'-a),  r^ev'. 


= 0,  rV  V', 


(98) 


where  n is  a constant  representing  the  number  of  fish  per  unit 
volume  and  V1  is  the  volume  of  the  sphere  of  radius  A.  With 
this  choice  for  n(?' ,a'),  Eq.  97  can  be  integrated  to  yield 
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G(r 1 , &')  = 


na 


y 


?'  eV 1 


(99) 


U) 


uQ 


= 0,  f’fV'. 


Now  the  Fourier  transform  of  the  average  scattered 
signal ^'|s(r, cu)\  and  the  Fourier  transform  of  the  average  signal 
within  the  school  •(^(r,  w))>  both  satisfy  appropriate  forms  of  the 
Helmholtz  equation,  Eq.  72. 

Since  G(?',ui)  = 0 exterior  to  the  school,  we  have 
y2<*s(r»>  + K2  <fS(t>)>  = 0,  (100) 

where  is  given  by  Eq.  94. 


For  the  wave  in  the  interior  of  the  school  we  have 


y2<(/(?»>  + k2(iu)  <*f(r,u>)>  = 0 , 


where 


k2(w) 


at2 

co 


4 Tina 


(101) 


(102) 


2 

To  obtain  this  expression  for  k (u>)  we  have  used  Eq.  73  and 


Eq.  99  for  r'eV'. 

Letting 


k = Y + ip, 


(103) 


where  Y and  {3  are  real,  we  find  from  Eq.  102  that 

2 4TTnaa^  (ua^-ou^) 

..  _ u.>  .•  v r ' 

Y 2 + "71  272  , 'TT’/ZTi — » 

co  (u,r"’tu  ) +U)  U)r/Q 


(104) 


4TTnau)  u)  /Q 

p = L 

, 2 2.2  , ,2  2/r.2 

(u>r  - to  ) +UJ  mr/Q 

Then 


(105) 


- -j  (y+\/y  +e2  ) 


(106) 


i r~2  21 
a - -y  (~V+\IY  -w  ) 


(107) 


where  K and  a have  been  defined  in  Eq.  22.  Thus,  the  wave 
propagating  through  the  school  experiences  absorption,  which  is 
frequency  dependent  since  a is  a function  of  uu,  and  dispersion 
since  k is  a function  of  uu.  The  phase  velocity  c of  the  wave 

r 

is  then  given  by  Eq.  24  and  in  general  this  will  differ  from  c . 
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We  note  for  future  reference  that  the  radial 
components  of  the  average  particle  velocities  <(u  (r,t»)),  and 
<(u^(r,  u>))>  for  the  scattered  and  internal  waves  respectively, 
are  given  in  terms  of  the  average  pressure  fields  by  the 
relations^ 


<us  (?,«;)>  - — i -§=  < 4s  (?,«.)>  , (108) 

1Poa, 

<uf (?,«»)>  - -i-  -§—  <(f(r»>  . ■ (109) 

ipfo) 

Likewise  for  the  incident  wave  we  have 

ui(?>u,)=JL_  _i_  f(r>)  . (110) 

iPouj  Br 

The  boundary  value  problem  then  consists  of 
finding  solutions  of  the  Helmholtz  Eqs.  100  and  101  that  satisfy 
the  boundary  conditions  (I)  at  the  surface  of  the  sphere  (r=A), 
the  acoustic  pressure  and  (II)  the  radial  components  of  the 
particle  velocity  must  be  continuous.  We  now  proceed  with  the 
solution  to  this  boundary  value  problem. 

27 

The  required  separated  solution  in  spherical 
coordinates  of  Eq . 100  is 

hv(K0r>  Vcos  0)  » v=  C111) 
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28 

where  Pv(cos  G)  are  the  Legendre  polynomials  and  hv(KQr)  is 
a spherical  Hankel  function  of  the  first  kind  defined  by 

hv(Kor)  = VKor)  + 1 W>  > (112) 

29 

j^(KQr)  and  nv(KQr)  being  a spherical  Bessel  function  and  a 

29 

spherical  Neumann  function  respectively. 

The  required  separated  solution  in  spherical 
coordinates  of  Eq.  101  is 

jv(kr)  Pv(cos  G)  , v=o,l,2 (113) 

Due  to  symmetry  these  solutions  (Eqs.  Ill  and  113)  are  independent 
of  cp. 

The  expansion  of  the  incident  wave  (,:L(r,  6 , uj) , given 
by  Eq.  96,  in  terms  of  the  eigenfunctions  of  the  Helmholtz 
equation  yields  x. 

. . iK  r cos  8 

r(r,6,(u)  = $L(m)  e 0 


v=o 


iv(2vfl)Pv(cos  0)jv(KQr)  . 


(114) 


The  scattered  wave  is  expanded  in  a similar  series 
with  coefficients  which  are  to  be  determined; 
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<;s(r,e,«.)>  -I^ib)  " iv(2vfl)Bv  Pv(ct)S  0)  hv(>.0r)  . (115) 

v=o 

Similarly,  for  the  internal  wave,  we  write 

CD 

< rtr,e,»»  = *V)  Y Avpv(e°s  6)jv(kr)  . (116) 

v=o 

The  coefficients  Av  and  Bv  are  found  by  applying 
the  boundary  conditions  I and  II  above.  The  first  boundary 
condition  requires  that 

^(A,e,to)+  <f(A,6,m)>  -<^(A,0,«j)>  , (117) 

and  the  second  boundary  condition  requires  that 

u^AjGjiu)  + ^us(A,e,cu)')  - </uf(A,6,w))>  • (118) 

Substituting  Eqs.  114,  115  and  116  into  Eq . 117 

yields 

jv(«0A)  + Bvhv(KQA)  = v(kA)  . (119) 

Similarly,  substituting  Eqs.  114,  115  and  116 
into  Eq.  118  after  having  made  use  of  Eqs.  108,  109,  and  110 


1 

1 
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yield 


jv<KoA)  + Bvhv<KoA)  = Z(^  > 


(120) 


where  the  prime  denotes  differentiation  with  respect  to  the 
argument  and  Z(w)  is  the  ratio  of  the  acoustic  impedance  of  water 
ZQ  to  that  of  the  fish  school  , i.e., 

Z 


Z(aO  = 


Zf((u) 


(121) 


where 


Zo  poco  » 


(122) 


Zf(O))  = , 


(123) 


Cf 


u)/k  . 


(124) 


Note  that  since  k is  complex,  the  speed  of  sound 
in  the  fish  school  c^  is  complex  and  consequently  the  ratio  of 
impedances  Z(iu)  is  complex. 


obtain 


Solving  Eqs.  119  and  120  simultaneously,  we 


jv(kA)j;(K0A)  - Z(»)  Jv(«„A)j;(kA)  1 
jv(kA)h^(KoA)  - Z(®)j;(kA)hv(K0A)  j ' 


(125) 
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Next  let  us  eliminate  the  derivatives  of  the 
Bessel  functions  in  Eq.  125  by  use  of  the  well  known  relation 


- -f  £v<5>  ■ £ \H-1 f " } ’ 


(126) 


where  f^(?)  can  or  * 


Using  this  result,  we  can  write  the  expression 


for  Bv  in  the  form 


N 


Bv  Mv  • 


(127) 


where 


Nv=  [vCK^-^l-hM^K^-j^C^A)]  jv(kA) 


+ Z(uj)  jv(*0A)  jy+l(1<A)  » 


(128) 


= [v0toA)'1(1-h)  hv(KoA)  " hvfl(KoA)]  jv(kA) 


+ Z(w)  hv(KQA)  3v+  i(kA)  , 


(129) 


and  we  have  put 


h = Pc/P£  . 


(130) 


Since  we  want  to  consider  the  scattered  wave  only 
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at  large  distances  from  the  fish  school,  we  can  use  the 

following  asymptotic  expansion  for  the  Hankel  function 

itc  r 

-v 

hv(KQr)  » , (Kor— >»)  . (131) 

iKor 

Thus  the  expression  (Eq.  115)  for  the  scattered 
wave  becomes  in  the  asymptotic  region 

iic  r 

(f  <r,8,u>)>  = «i(w)S(e,a))  ^ , (132) 


where 

S(e,cu)  = (i^0)-1  y (2v+l)Bv(w)Pv(cos  6). (133) 

v=o 

Consequently,  the  expression  for  the  scattered 


signal  becomes 

. ✓ 

« -iu)(t-r/co) 

^Ys(r,G,t))>  = Re  J ^(ui)  S(9,m)  - dm  . (134) 

o r 

Now  if  we  write 

S(9,m)  = Sj  + i S2  , (135) 


where  S-^  and  S2  are  real,  Eq.  134  becomes 
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<V(r,e,t))>  - rjrp-  J*  ^ 1(a;)  [s.^8 , u))cos  u>t*+S2(9 , u>)sin  cut*]  dx,  (136) 


where  t*  is  the  retarded  time 


t*  - t - r/cQ  • 


(137) 


Finally,  let  us  calculate  the  scattered  signal 
in  the  Born  approximation.  The  Fourier  transform  of  the  signal 


is  given  by  Eq.  76 


iK  r-r 


As(r,e,u»)>  = J G(?»^(?',«»)  ~ — — 

N r-r 


(138) 


Let  (r,9,cp)  be  the  spherical  coordinates  of  r 
and  (r'jG'jCp')  those  of  r ' . Further,  let  © be  the  angle  between 
r and  r'  (c£.  Fig.  3). 


fl \ ' 
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Fig.  3 
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Using  Eq.  99  for  G(r',ui)  and  Eq . 96  for  ^(r ' , u;) , 
the  expression  for  the  scattered  wave  given  by  Eq.  138  becomes 


TT 


2tt 


<^s(r,0,u)))  = G(uj)^1(aO  J dr'  J d©  * J dcp'  r'^  sin  8' 


iK  r ' cos  8'  ^'*'Ko ! r r 


x e 


r-r ' 


(139) 


->i 


r-r 


-[ 


Now 


r + r'2 


2rr'  cos  ©j 


1/2 


(140) 


where 


cos  @ = cos  6 cos  0'  + sin  8 sin  0'  cos  (cp-cp1)  . 


(141) 


Since  we  want  to  consider  the  scattered  wave 
only  in  the  far  field,  we  have  r'«  r.  Consequently,  we  can 
expand  Eq.  140  retaining  only  first  order  terms: 

r r'2  r'  ]l/2 

^ 1 + - 2 ~ cos  Bj 


?-?'!  = r 


•[i- 


cos  © 


r-r'  cos  © 


(142) 


While  it  is  necessary  to  retain  the  first  order 
term  r'  cos  © as  far  as  the  factor  |r-r'|  which  appears  in  the 
phase  is  concerned,  it  is  not  necessary  for  the  factor 
|?-r  'j  which  appears  in  the  amplitude.  Consequently  for  the 
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amplitude  we  can  write  |r-r'|«  r. 

Moreover,  we  wish  only  to  consider  the  back- 
scattered  wave  in  the  direction  6 = TT,  so  that  cos  © = - cos  6 1 
from  Eq.  141.  Collecting  all  of  these  results  allows  us  to 
write  Eq.  139  as 

0 . r lKorlA  TT  2TT  o 2iK  r 'cos  G' 

/tyS(r,  TT,  u>)^>  = G(  u.)  § ^ ( uj)  ~ — J dr'  j'  dO  1 J'  dcp'r1  sin  0'  e 
^ L ! o o 


(f(r,n,w))  = G(iw)*r(«))^ 


<^s(r,TT,u>))>  = 2lT  G(w)ti(iv) 


iK  r . „ 2iK  r'  cosG' 

O A II  o O 

— J dr'  j dG'r'Z  sin  6'  e 

o o 


<f  (r,TT,«>)>  = 2T TG(«J)  ^(iv) 


iKor  A 

- J sin  (2k  r')  r'dr' 

K r o ° 

o 


<'>S(r,TT,aj))  = Jig ( !Ll) [ s in  ( KqD) - KpD  cos(KoD)] 

2k3 

o 


iK  r 

i ° (143) 


where  D = 2A. 


The  scattered  pulse  in  the  Born  approximation 


then  becomes 


i 05  -i  wt 

{Ys(r,TT,t))  =—  Re  J ^('S(r,TT,u)))  e dw 


<^s(r,TT,  t )^>  - Re  J 


J [sin(KoD)-KoD  cos(KoD)]  e ^ 

o k3 
o 


daj  (144) 


where,  as  before,  t*  is  the  retarded  time  defined  by  Eq.  137. 


G('.u)  - + i G2  , 


(145) 


where  from  Eq.  99,  we  find  that 


2/2  2S 

naa;  (u-  -to  ) 

G1  = /~2  272“ T 2 ,„2  ’ 

( our  - tu  ) + tt>rti  /Q 


(146) 


g2  = 


O 

nato  u)r/Q 

i 2 2 2 /n2 

(u)r-d  ) + tord  /Q 


(147) 


then  we  can  write  Eq.  144  as 


00 

<V(r,TI,t;*)>  = 4?  J"  [sin(E0D>-KoD  cos(KoD>] 

O K 


x ^ G-^  cos  tit*  + G2  sin  tot*]  dio 


(148) 


This  is  the  average  of  the  scattered  signal  in 
the  Bom  approximation.  Note  that  the  scattered  signal  is 
directly  proportional  to  the  average  number  of  scatterers  per 
unit  volume  n as  can  be  seen  from  Eqs  146  and  147. 


I 
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D.  SCATTERING  FROM  A LAYER  OF  FISH  WITH  PLANE  BOUNDARIES 

Consider  a random  collection  of  fish,  all  with 
identical  scattering  properties  and  uniformly  distributed 
within  a volume  bounded  by  two  parallel,  infinite  planes.  Let 
us  choose  a coordinate  system  so  that  the  planes  are  parallel  to 
the  xy-plane  and  let  one  plane  intersect  the  z-axis  at  z = h 
and  the  other  intersect  at  z = h + Ah  (Ah  > 0) . The  geometry 
is  illustrated  in  Fig.  4. 
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In  reference  to  Fig.  4,  media  I and  3 are  water  while 
medium  2 is  a mixture  of  water  and  fish.  Since  the  same 
as sump tions  will  be  made  for  the  layer  of  fish  as  were  made  for 
the  spherical  school,  full  use  will  be  made  of  the  results  in 
Section  C without  restating  or  rederiving  them  here. 

Let  be  the  density  and  c^  the  sound  velocity  in 
the  i medium.  Then,  from  what  was  said  above,  = P3  = PQ> 

— , c ^ — * C3  ~ ~ cQ  and  . 

Let  Y^  (r , t)  be  a plane  pulse,  traveling  in  medium  1, 
which  is  incident  on  the  layer.  Since  the  boundary  of  the  layer 
is  a plane,  it  is  more  convenient  to  use  the  Fourier  plane  wave 
expansion  for  Y^  (?, t)  than  the  general  expansion  given  by 
Eq.  19.  To  obtain  the  plane  wave  expansion  from  the  general 
expansion  (Eq.  19),  one  simply  lets  the  Fourier  transform 
^ (?,«>)  be  a plane  wave: 


(?,«•>)  = 51(m)  e 


i£i 


(149) 


where  $x(w)  is  the  amplitude  of  the  wave  and  depends  only  on 
the  frequency  w.  For  the  case  of  a CW-pulse,  41(o>)  is  given 
by  Eq.  95. 

The  vector  is  called  the  wave  vector.  Its  direction 
is  that  of  the  normal  to  the  wavefront  and  its  magnitude, 
denoted  k-^,  is  given  by 


k,  1" 

1 0 c0  *0 


(150) 


Let  n be  the  unit  normal  vector  to  the  wavefront  and 

let  n , n , n be  the  components  of  n along  the  x,  y and 
x y 2 

z-axes,  respectively.  Since  n is  a unit  vector,  the  components 
must  satisfy  the  relation 


2 2 2 , 
nx  + "y  + nZ  = 1 ' 


(151) 


Now  we  can  write  in  the  form 


kl  = kln  » 


(152) 


and  the  x,  y and  z components  of  k^,  denoted  k].x>  ^lz’ 


in  the  form 


klx  - klnx>  kly  - klny  klz  ‘ klnz  > 


(152) 


so  that 


.2  . , 2 , , 2 ,2,2.  2 , 2N  ,2 
klx+  kly  + kU  - kl  <"x  + ny  + nz>  ' kl  • 


(153) 


Now  consider  the  situation  illustrated  in  Fig.  4 where 

the  unit  normal  n to  the  incident  wavefront  lies  in  the  xz-plane 

and  makes  an  angle  9p  with  the  positive  z-axis.  In  this  case 

n = sin  6,  , n = 0.  n = cos  8-.,  and  Eq.149  becomes 
x 1*  y ’ z l’ 


^ (?,u.')  = *i(uj)  e 


^lx*  + klyy  + klz  X) 


- .!(.)  elkl(V  + 

. ik,  (x  sin  6,  + z cos  0.) 
= «i(u:)  e 1 1 


(154) 
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Thus  the  Fourier  integral  representation  (Eq.  19) 


for  Y^(r, t)  becomes 


Y^(x,z,t)  = ^ Re  / * (uj)  e 


. ik,  (x  sin  0,  + z cos  0,)  -ii 

i / \ _ l i.  J. 


duj  . 

(155) 


A similar  Fourier  plane  wave  expansion  can  be  made 
for  the  average  of  the  signal  scattered  back  into  medium  1: 


<jY®(x,z,t))  = yf  Re  J'  R( iu)e 


ik^(x  sin  0^  - z cos  0^)  -irnt 


dm  , 


(156) 


where  R(m)  is  the  reflection  coefficient  to  be  determined  by 
boundary  conditions. 

Let  ^Y^(x,z,t)]>  denote  the  total  sound  field  in 
medium  1 so  that 


^Y1(x,z,t))  = Y^(x,z,t)  + <Y^(x,z,t))  . 


(157) 


We  can  represent  the  average  of  the  total  sound  field 
in  medium  2,  ^Y2(x,z,t)^  , and  in  medium  3,  <^Y3(x,z,t))>  , in  a 


similar  manner: 


CO 

^Y2(x,z,t))  = \ Re  J A2(«0  e 


ik2(x  sin  @2  + z cos  02)-imt 


. to  ik9(x  sin  09  - z cos  8o)-iu.’t 

+ Re  J B2(uj)  e dm  , 


(158) 
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•,  00  iko(x  sin  G-  + z cos  Go)-iait 

^Y3(x,z,  t))  = jj  Re  J T(u)>  e J J J 


<Y3(x,  z,  t))  = Re  J T(o)>  e J J duo. 

0 (159) 

Here  k2  = k , k3  = Kq  and  02  Is  the  angle  of  refrac- 
tion in  medium  2 while  6^  is  the  angle  of  refraction  in  medium  3 
(Cf.  Fig.  4).  The  amplitudes  A2(m),  B2(u')  and  T(co)  of  the 
respective  waves  are  to  be  determined  by  the  boundary  conditions, 
The  first  integral  in  Eq.  152  represents  the  total  sound  field 
which  is  propagating  in  the  direction  with  k2Z  > 0 and  the 
second  integral  represents  the  total  sound  field  which  is 
propagating  in  the  direction  with  ^2z  < 

Since  the  normal  to  the  boundaries  is  parallel  to 
the  z-axis,  the  normal  components  of  the  particle  velocities 


are  given  by 


/u  (x  z «))  = -T-i — ~ /ty  (x  z , uu))> 
\ mv  ’ * ''  ip  ® dz  \ mv  * * '/ 


(160) 


for  m = 1,2,3. 

The  boundary  conditions  are  then  that  the  pressure 
and  the  normal  components  of  the  particle  velocities  must  be 
continuous  across  both  planes  comprising  the  boundary. 


Continuity  of  pressure  at  z = h: 


<*!_  = <^2(x,h,u)))>  . 


(161) 


Continuity  of  pressure  at  z = h + Ah: 


<^^2 (x > h + A*1*  wY}  ~ ^i3(x,h  + Ah , m) y . 


(162) 


Continuity  of  the  normal  component  of  particle  velocity  at 
z = h: 


^(x^h.uj))  =<u2(x,h,a >)^>  . 


(163) 


Continuity  of  the  normal  component  of  particle  velocity  at 


z = h + Ah; 


<u  2 (x > + Ah,  u>))>  =(u  j(x,h  + Ah,  m))>  . 

Using  Eqs.  155-158,  the  first  boundary  condition 


(164) 


(Eq.  161)  becomes 


ik,  (x  sin  0,  + h cos  0,)  ik,  (x  sin  8,  - h cos  6,) 

>V')  e 1 1 1 + R(u)  e 1 1 1 


ik0(x  sin  09  + h cos  0?)  ik9(x  sin  0?  - h cos  ©2) 

A,(w)  e 2 2 2 + B2(u.)  e 2 

(165) 


Rearranging  the  terms  in  this  equation,  allows  us  to 


write  it  as 


ix(k,  sin  - k£  sin  02)  r ik2hcos02  ik2hcos02' 

> = A2e  + B2e  J 


' . ik,h  cos  0.  -ik,h  cos  0,  -1 

X }le  1 1 + Re  1 


(166) 


Since  the  right  side  of  this  last  equation  is 
independent  of  x,  the  left  side  must  also  be  independent  of 
x.  Consequently,  we  must  have 
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This  is  just  Snell's  law  of  refraction  for  media  1 and  2. 

Using  Eqs.  158  and  159,  the  second  boundary  condition 
(Eq.  162)  becomes 


Combining  Eqs.  167  and  170,  yield 

sin  0^  = k2  sin  ^ = k^  sin  0^  (171) 
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Using  Eqs . 155-158  and  160,  the  third  boundary 


condition  (Eq.  163)  becomes 


ik,  ( x sin  6^  + h cos  0-^) 

_i [ ik,  cos (0-,  )SL  e 

TIEpJ  [ 1 l 

ik,  ( x sin  G1  - h cos  0-^ 


- ikx  cos^^R  e 


_ _1_  ik0  cos  (6o)A2 

rujp2  L 2 L 


ik9( x sin  S2  + ^ COS  ^2^ 


ik0( x sin  02  - h cos  02) 


- ik2  cos(02)B2  e J 

Finally,  using  Eqs.  158-160,  the  fourth  boundary 


(172) 


condition  (eq.  164)  becomes 


lK-2(,  * u2 

1 rik2  cos (@2^2  e 


( x sin  G2  + h cos  02+Ah  cos  02) 


imP2  l 


i k2  cos(02^®2  e 


ik2(xsin  Q2  " h cos  02«Ah  cos  S2)l 


i—  iko  cos(S~  )T  e 
iwP3  I 3 3 


ikj  ( x sin  03+h  cos  G3+£h  cos  ®3) 


.(173) 
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If  we  use  Eq.  171  and  the  fact  that  kj-rc/c^  Eqs. 
165,  168,  172  and  173  can  be  simplified  to 

ikh  cos  e -ikth  cos  ik2h  cos  e 'lk2h  COs62 

1 +Re  -A9e  2 


,(174) 


ik0(h+Ah)cos  89  -ik2(h+Ah)cos  82 


+B2e 


iko (h  cos  03+  Ah  cos  83) 


= Te 


(175) 


$i  ik^h  cos  0^ 


-ikjh  cos  6L  _ A2  ik2h  cos  82 

It e = z7e 


B9  "ik2h  COS  82 

^ A 


(176) 


A,  ik2(h«h)cos  92  _ _B2_  -ik2(h«h)cos  02 


ik3(h+Ah)  cos  83 


where  we  have  defined  the  impedance  Z±  as 


(177) 


_ pici 
zi  “ cos 


(178) 
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Eqs.  175-178  represent  a system  of  four  equations 
for  the  four  unknown  amplitude  functions  R,  A2>  B2  and  T. 
When  they  are  solved  simultaneously  for  R they  yield 


R(w)  = 


V21  + V32  e 


2iVh  cos  ®2  " 21k,h  cos 
2ik0Ah  cos  OT 


1 + V21V32  e 


(179) 


where  we  have  put 


Z.  + Z. 
1 1 


(180) 


Since  we  are  only  interested  in  the  reflected  wave, 
the  other  amplitudes  will  not  be  determined. 

The  result,  Eq.  179,  can  be  put  into  a more 
convenient  form  for  computation.  If  we  let 


D.  . 
ij 


2 Z. 

1 ■■  J 


Zi  + Zj 


(181) 


and  expand  the  denominator  of  Eq.  179  in  a power  series,  we 
obtain 


-2ik,h  cos  ei  2ik ,4b  cos  6, 

( ,i,  , “ V21  + D12  V32  D21  e 

5 (uc) 


+d12  V32  V12  V32  D21  e 


4ik2Ah  cos  ©2 


+.  ..,  (182) 


or  written  more  compactly, 


. -2ik,h  cos  G,  0.,  0 

R(m)  e 1 1 . _2rk2Ah  cos  G2 

V2i+Di2  V32  D21  e 


^(w) 


x / V12V32  e 


2ik2Ah  cos  02  i n 


(183) 


This  last  form  (Eq.  183)  for  the  reflection  coefficient 

31 

has  an  interesting  physical  interpretation.  It  is  the  resolution 

of  the  total  reflected  wave  into  individual  waves  which 

represent  the  various  multiple  orders  of  reflection  at  the 

boundaries  of  the  layer.  The  coefficient  V..  can  be  interpreted 

as  the  reflection  coefficient  of  the  boundary  between  the  l 

and  j media  for  a wave  that  is  traveling  from  the  j ^ 
th 

medium  to  the  i medium.  Likewise,  the  coefficient  D..  can 

be  interpreted  as  the  transmission  coefficient  of  the  boundary 
tli  th 

between  the  l and  j media  for  a wave  that  is  traveling  from 

th  th 

the  j medium  to  the  ic  medium. 
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Using  this  interpretation  of  and  D^.  , let  us 

consider  the  individual  terms  in  Eq.  182.  The  first  term 

V21  is  ^e  complex  amplitude  of  a wave  that  has  been  reflected 

once  from  the  boundary  between  media  1 and  2.  The  second  term 
2ik2Ah  cos  62 

D21  V^2  e is  t^ie  comP^-ex  amplitude  of  the 

wave  which  has  been  transmitted  from  medium  1 to  medium  2 
(D21) > passed  through  medium  2,  was  reflected  at  the  boundary 
between  media  2 and  3 (V^)}  passed  back  through  the  layer  and 
back  into  medium  1 again  (D-^)*  The  additional  phase  shift 
which  occurs  in  the  exponential  is  simply  due  to  the  slant 
distance  across  the  layer  which  is  just  (Ah  cos  0 2) • The 
succeeding  terms  may  be  interpreted  similarly. 
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